Ab initio study of the giant ferroelectric distortion and pressure induced spin-state 

transition in BiCoOa 
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Using configuration-state-constrained electronic structure calculations based on the generalized 
gradient approximation plus Hubbard U method, we sought the origin of the giant tetragonal ferro- 
electric distortion in the ambient phase of the potentially multiferroic material BiCoOa and identified 
the nature of the pressure induced spin-state transition. Our results show that a strong Bi-0 co- 
valency drives the giant ferroelectric distortion, which is further stabilized by an xy-type orbital 
ordering of the high-spin (HS) Co'^"'" ions. For the orthorhombic phase under 5.8 GPa, we find that 
a mixed HS and low-spin (LS) state is more stable than both LS and intermediate-spin (IS) states, 
and that the former well accounts for the available experimental results. Thus, we identify that the 
pressure induced spin-state transition is via a mixed HS-(-LS state, and we predict that the HS-to-LS 
transition would be complete upon a large volume decrease of about 20%. 
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PACS numbers: 75.30.-m, 71.20.-b, 71.27.-|-a, 71.15.Mb 



I. INTRODUCTION 

Multiferroic materials, having coexisting magnetism 
and ferroelectricity, are of great technological and funda- 
mental importance, 1"'^ given the prospect of controlling 
charges by applying magnetic fields and spins by voltages. 
BiCoOa was recently synthesized by a high-pressure (HP) 
technique,^ and it has been suggested to be a promising 
multiferroic material by Uratani et al.^ and by Ravindran 
et al^ both through first-principles Berry-phase calcula- 
tions. BiCoOa has a giant tetragonal lattice distortion 
of c/a = 1.27 with remarkable off-center atomic displace- 
ments (see the inset of Fig. 1), and it is an insulator 
having C-type antiferromagnetism below 470 K — the an- 
tiferromagnetic (AF) ab layers stacking ferromagnetically 
along the c axis;^ 

It was proposed^i^ that the giant tetragonal distor- 
tion originates from lifting of the orbital degeneracy of 
the high-spin (HS, 5'=2) Co'^^ ions and is stabilized 
by the subsequent xy-type ferro-orbital ordering. Note 
that orbitally degenerate transition-metal oxides quite 
often display an orbital ordering (00) but ferroelectric 
(FE) materials out of them are rare, as ferroelectricity 
and magnetism seem, and actually in most cases, to ex- 
clude each other i^iiS Therefore, the proposed mechanism 
for the giant FE distortion appears not straightforward. 
Using fixed-spin-moment density-functional calculations, 
Ravindran et al.^ predicted that there is a giant magneto- 
electric coupling in BiCoOa: an external electric field (or 
a small volume compression of ~5%) can induce a strong 
magnetic response by changing the magnetic Co^^-HS 
state in the FE phase into a nonmagnetic low-spin (LS, 
5=0) state in a paraelectric (FE) phase. A corresponding 
HS-insulator/LS- metal transition was also suggestedi^iii 
In sharp contrast, a very recent HP study^ showed that 



BiCoOa even under 6 GPa with a large volume decrease 
of 18% is still semiconducting. Note however that con- 
troversial spin states, both LS and intermediate-spin (IS, 
5=1), were suggested for the HP phase.— 

In the present work, we seek the origin of the giant 
tetragonal FE distortion in the ambient phase of BiCoOa 
and identify the nature of the pressure induced spin- 
state transition, using two sets of configuration-state- 
constrained GGA-I-J7 (generalized gradient approxima- 
tion plus Hubbard U) calculations. Our results show 
that the giant tetragonal distortion is driven by a strong 
Bi-0 covalency (rather than by the aforementioned lift- 
ing of orbital degeneracy) and is further stabilized by an 
xy-type 00 of the HS Co'^^ ions. Moreover, we find that 
the pressure-induced spin-state transition is via a mixed 
HS-I-LS state, which accounts for the available experi- 
mental results consistently and disproves a recent pre- 
diction of a readily switchable HS-LS transition. 



II. COMPUTATIONAL DETAILS 

We used the structural data of BiCoOa measured by 
the neutron power diffractiouji Our calculations were 
performed using the full-potential augmented plane- wave 
plus local-orbital code WIEN2k.^^ The muffin-tin sphere 
radii were chosen to be 2.3, 1.9, and 1.4 bohr for Bi, Co, 
and O atoms, respectively (1.0 bohr for O when calculat- 
ing the E-V curves shown in Fig. 6). The cutoff energy 
of 16 Ryd was set for the plane- wave expansion of in- 
terstitial wave functions, and 1200 k points in the first 
Brillouin zone for the ambient structure with one formula 
unit (f.u.) and 300 k points for the HP phase with 4 f.u. 

Plain GGA [or local-spin-density approximation 
(LSDA)] calculations^i^iiiiH seem to qualitatively re- 
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FIG. 1: (Color online) E vs c/a curve calculated by GGA+!7 
for the Co'^'''-LS relaxed structures of BiCoOs with c/a = 1- 
1.3 (in a step of 0.05) including 1.27 (expt.). The inset shows 
a unit cell of the experimental tetragonal structure. 



produce the C-type AF and insulating ground state of 
BiCoOa in the ambient phase, which was ascribed to the 
strong Hund-exchange stabilized HS state of the 00-^+ 
ions (and thus the AF order and narrow bands) and to 
the well split-off xy-singlet orbital. However, the band 
gap of 0.6 eV and the Co'^+ spin moment of 2.4 fiB given 
by the GGA/LSDA calculations are much smaller than 
the experimental values of 1.7 eV and 3.2 /is.'*'^^ More- 
over, a recent prediction of an insulator-metal transition 
in BiCoOa upon a volume decrease of ~5% made by 
GGA/LSDA calculations^iii has already been disproved 
by a very recent HP study, '' which shows that BiCoOa is 
still semiconducting even under 6 GPa with a large vol- 
ume decrease of 18%. As seen below, the experimental 
values of both the band gap and the Co'^"'" spin moment 
are well reproduced by our GGA+U calculations. Note 
also that BiCoOa has an apparent AF order up to 470 Ki^ 
All these suggest that BiCoOs should rather be catego- 
rized as a Mott insulator, with its band gap determined 
primarily by Hubbard U, i.e., strong correlation of the 
Co 3c? electrons. 

To account for the strong electronic correlation^i^ we 
have carried out GGA-I-J7 calculations throughout this 
paper. In particular, we used the configuration-state- 
constrained GGA+U method, ^^"^'^ which allows us to ac- 
cess different spin and orbital configuration states of the 
concern by initializing their corresponding density ma- 
trix and then doing self-consistent electronic relaxation. 
This method is quite useful for study of the spin and or- 
bital physics present in transition-metal oxidesj^^"— All 
the results shown below are obtained with U=Q eV and 
Hund exchange J— 0.9 eV. We note that our test calcu- 
lations using other reasonable U valueSfi^iiSiiiS U=5 and 
7 eV, gave qualitatively the same results. 



TABLE I: The total energies of BiCoOa (in unit of eV/f.u.) 
relative to the hypothetical cubic-structure LS state calcu- 
lated by GGA+U : the experimental tetragonal phase assum- 
ing an LS state with atomic relaxations (the third row), as- 
suming an HS state in the LS relaxed structure (the fourth 
row), and assuming the HS state with further atomic re- 
laxations (the fifth (last) row). The optimized atomic z- 
coordinates are also shown. Note that the structural data of 
the atomically relaxed HS state are in good agreement with 
the experiment^ (see the inset of Fig. 1). See more discussion 
in the main text. 
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-1.47 






Tetragonal HS-relaxed 


-1.97 


0.5631, 


0.1912, 0.7237 



III. RESULTS AND DISCUSSION 

We first seek the origin of the giant tetragonal distor- 
tion by starting with our calculations assuming an ideal 
cubic structure (space group PmSm) having the same 
volume as the experimental ambient tetragonal structure. 
Then we changed the c/a ratio of the tetragonal struc- 
ture (keeping the volume unchanged) in our calculations. 
All those calculations were carried out by setting a hy- 
pothetical LS state of the Co^+ ions, and by doing a full 
electronic and atomic relaxation for the cases oi c/a — 
1.05-1.3 in a step of 0.05, and for the experimental 1.27 
as well. All the solutions are insulating. As the LS Co'^"'" 
has a closed sub-shell ^ (thus no orbital degeneracy) and 
is an isotropic ion, we can use this set of LS-constrained 
GGA-I-J7 calculations as a computer experiment to probe 
mainly the Bi-0 covalent effect. 

We show in Fig. 1 the calculated total energies as func- 
tion of the c/a ratios, and one can immediately find that 
the c/a ratio in the LS equilibrium state is close to 1.2 
(about 1.18), already indicating a large tetragonal dis- 
tortion. The corresponding energy gain, relative to the 
hypothetical cubic structure (c/a = 1), is 0.96 eV/f.u. 
For the LS relaxed structure with the experimental c/a 
= 1.27, the energy gain is 0.90 eV/f.u. We list in Table I 
(see the third row) the optimized atomic z-coordinates of 
the LS state with the experimental c/a ratio, and we find 
that the experimental C0O5 coordination is well achieved 
[as seen from the optimized Co-0 bondlengths, Co-01: 
1.893 Ax 1, 2.832 Ax 1 (much larger) ; Co-02: 1.919 Ax4] 
even in the presence of the Co^+ LS state without or- 
bital degeneracy. This is also the case for c/a — 1.2 with 
another optimized Co-0 bondlengths, 1.881 Axl, 2.677 
Axl, and 1.943 Ax4. We plot in Fig. 2(a) a charge den- 
sity contour of the hypothetical cubic structure in the 
(100) plane, and in Fig. 2(b) that of the LS relaxed 
state with the experimental c/a ratio. Fig. 2(a) shows 
nearly spherical charge densities around both the Bi and 
O ions, indicative of Bi^+-0^~ ionic bonds in the hypo- 
thetical cubic structure. In contrast. Fig. 2(b) clearly 
shows directional Bi-0 covalent bonds, which are also in- 
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dicated the orbitally resolved density of states (see Fig. 
3 and the discussion below). All these suggest that the 
Bi-0 covalency drives the giant tetragonal FE distortion 
and causes remarkable off-center atomic displacements, 
no matter that the Co^+ ions are in this hypothetical LS 
state or in the real HS state (see below) . This is similar to 
the stereochemical mechanism of the Bi'^"'' 6s^ lone-pairs 
proposed for the highly distorted perovskite manganite 
BiMnOa^ 



TABLE II: The relative total energies AE (meV/f.u.), Co^+ 
spin moments M {hb), and band gap Eg (eV) of BiCoOa in 
the 5.8 GPa phase calculated by GGA+C/ for the LS-relaxed, 
IS-relajced, and (HS+LS)-relaxed structures. 
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FIG. 2: (Color online) Charge density contour (0.01-0.1 e/A^ 
in a step of 0.01 e/A^) of BiCoOa in the (100) plane of (a) the 
hypothetical ideal cubic structure and of (b) the LS-relaxed 
structure with the experimental c/a ratio. The directional 
Bi-0 covalency is apparent in (b). 
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FIG. 3: Partial density of states of BiCoOa in the C-type AF 
ground state of the ambient phase. 

Furthermore, we study the effect of the xy-type 00 
of the HS Co'^"'" ions on the giant tetragonal distortion. 
Owing to the effective C0O5 pyramidal coordination al- 
ready present in the LS relaxed tetragonal structures as 
described above, the crystal field (particularly the t2g-eg 
splitting) is weak and thus less important than the Hund 
exchangej2i As a result, an HS state with an 00 of the 
well split-off xy orbital is the ground state and more sta- 
ble than the hypothetical LS state by 0.57 eV/f.u. for the 



experimental c/a — 1.27 (see the third and fourth rows in 
Table I) and by 0.59 eV/f.u. for c/a — 1.2. In this sense, 
the xy 00 is a consequence of the Bi-0 covalency driven 
tetragonal distortion and the associated C0O5 pyrami- 
dal coordination. Moreover, our calculations doing struc- 
tural optimization show that an adjustment of the lattice 
to this xy 00 of the HS Co''+ ions changes the c/a ratio 
to 1.28 (as compared to 1.18 for the LS equilibrium state, 
see Fig. 1). This agrees very well with the experimental 
c/a — 1.27. Correspondingly, the experimental atomic 
parameters are also well reproduced, see the last row in 
Table I. Furthermore, we show in Fig. 3 the density of 
states of BiCoOa in the C-type AF ground state having 
the xy 00 of the HS 00"^+ ions. The calculated band gap 
of 1.98 eV and the Co'^^ spin moment of 3.01 /is are also 
in good agreement with the experimental values of 1.7 
eV and 3.24 ^b^^ Hybridizations between the Bi 6s6p 
and O 2p orbitals are also apparent and evidence again 
the Bi-0 covalency, although the magnitude of their re- 
spective density of states is much underestimated within 
the muffin-tin spheres as those orbitals are spatially quite 
spread. In a word, all above results allow us to conclude 
that the giant tetragonal FE distortion of BiCoOa origi- 
nates from the Bi-0 covalency (rather than from lifting 
of the orbital degeneracy of the HS Co^+ ions^'^^^) and 
is further stabilized by the subsequent xy-type 00. 

We now identify the nature of the pressure induced 
spin-state transition in BiCoOa, using another set of 
constrained GGA-fC/ calculations for the orthorhombic 
structure measured at 5.8 GPa, by assuming the LS, IS, 
and the mixed HS-I-LS states, respectively, and by doing 
a full electronic and atomic relaxation for each case. 

As seen in Table II, the mixed HS-I-LS state has the 
lowest total energy, and the LS (IS) state lies above it 
by 91 (158) meV/f.u. Those results suggest that either 
the pure LS or IS state present in the 5.8 GPa phased 
is not the case. Reversely, if the pure LS state were 
present in the 5.8 GPa phase, it would give rise to a 
change of the spin state from the pure HS state in the 
ambient phase, AS'=2, being in disagreement with the 
observed AS'=1.^ Moreover, an absence of the IS state is 
also not surprising, as (1) BiCoOa in the 5.8 GPa phase 
is free of a Jahn- Teller distortion (that is expected for 
the localized IS Co^'^); (2) a half-metallic band struc- 
ture of the IS state (see Fig. 4(a)) disagrees with the 
measured semiconducting behaviorfi and (3) up to now 
a definite example of the insulating IS state appears still 
lacking, and even in the layered perovskites LaSrCo04 
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FIG. 4: (Color online) The total and Co 3d density of states 
of BiCoOs in the 5.8 GPa phase calculated by GGA+U for 
the relaxed IS (a), LS (b), and mixed HS+LS (c) structures, 
respectively. 



and Lai.5Sro.5Co04 both having a strong tetragonal elon- 
gation of the Co^+Og octahedra, the IS state turns out 
not to be the ground state either^iS, despite an IS state 
might be intuitively expected. 

As such, the mixed HS+LS state is most probably 
present in the 5.8 GPa phase. An ideal 1:1 configura- 
tion of the mixed HS-I-LS state, with an average 5=1, 
well accounts for the observed change of the spin state*2 
Moreover, a G-type order of the HS and LS Co^+ ions 
(each HS Co ion is surrounded by six LS Co ions, and vice 
versa), due to a bigger /smaller size of the HS/LS Co^+ 
ions, could help to gain an elastic energy, and the resul- 
tant Co-0 bondlengths are calculated to be 1.980 (1.930), 
1.976 (1.890), and 1.948 (1.906) A for the HS (LS) Co^+ 
ions along the local xyz axes. Note, however, that a long- 
range G-type order of the HS and LS Co^+ ions is hard 
to establish, as only single transition-metal species in a 
single charge state and in the identical octahedral coor- 
dinations, i.e., solely Co'^"'' ions are involved. Thus, the 
average Co-0 bondlengths of the mixed and disordered 
HS-I-LS state (possibly with a short-range order due to 
a partial release of the lattice elasticity) are also in good 
agreement with the experiment. - Furthermore, the cal- 
culated insulating gap of BiCoOa is reduced from 1.98 eV 
in the C-type AF state of the ambient phase to 0.84 eV 
in the mixed HS-I-LS state (see Fig. 4(c)), which quali- 
tatively accounts for the decreasing resistivity of BiCoOs 
under pressure.^ 

Note that when considering thermal excitation of elec- 
trons into the initially empty conduction band and holes 
left in the valence band for a nominally stoichiometric 
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FIG. 5: Electron hopping (a) from an HS Co^^ ion to a neigh- 
boring HS Co^+ and (c) from an LS Co='+ to an LS Co*+, but 
a suppressed electron hopping (b) from an HS Co^^ to an LS 
Co*+ due to a spin blockade. 



material, these electron excitations would in a localized 
picture correspond to HS Co^"*" states and holes to LS 
Co4+. They behave like HS Co^+ and LS Co^+ "impu- 
rities" in the matrix of the Co'^"'" ions. As a result, in 
the ambient phase of BiCoOa having the HS Co'^"*" ma- 
trix, only the HS Co^+ "impurities" can transfer their 
minority-spin t2g electrons to the neighboring HS Co^"*", 
without changing the configuration states (the initial and 
final states are the same), see Fig. 5(a). However, a hole 
hopping from the LS Co^"*" to the HS Co^"*" is significantly 
suppressed (see Fig. 5(b)), due to a cost of the Hund 
exchange energy associated with a large change of the 
spin states which is referred as to a spin-blockade mach- 
anismi^ In contrast, in the HP phase having the mixed 
and disordered HS-I-LS Co'^"'" matrix, a charge hopping 
can take place both between the HS Co^"*" and HS Co'^^ 
(Fig. 5(a)) and between the LS Co^+ and LS Co^+ (Fig. 
5(c)). This could also account for the decreasing resis- 
tivity of BiCoOa under pressure. 

Our above results show that even in the 5.8 GPa phase 
of BiCoOs with a large volume decrease of 18%, the HS- 
to-LS transition is not yet complete, and the system is 
most probably in the mixed HS-I-LS insulating state, but 
not in the pure-LS metallic state which was predicted by 
the previous GGA/LSDA calculations for BiCoOs upon 
a volume decrease of >5%i^iii As we discussed above, 
BiCoOa is rather a Mott insulator, and thus the previ- 
ous GGA/LSDA results are somewhat questionable: par- 
ticularly the predicted metallic solution was already dis- 
proved by a very recent experimentji Recent fixed-spin- 
moment calculations even predicted that BiCoOs could 
have a giant magnetoelectric coupling with a readily 
switchable HS-LS transition associated with an electric 
field driven FE-PE transitioni^ We note that the predic- 
tion may simply be an artifact of the fixed-spin-moment 
calculations, as (1) LSDA or GGA (it was mentioned as 
a density- functional method in Ref®) is not well suitable 
to describe this Mott insulator; (2) most probably those 
calculations were carried out in a wrong ferromagnetic 
metallic state; and (3) their metallic solutions blurred the 
distinction between the different spin and orbital mul- 
tiplets of the concern and thus suppressed significantly 
their level splittings and particularly the HS/LS split- 
ting: the fixed-spin-moment calculations showed that for 
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FIG. 6: E-V curves calculated by GGA+f/ for the C-type 
AF state of BiCoOs in the ambient phase, and for the mixed 
HS+LS state and the LS state under high pressures. The 
lattice volume of the ambient phase is slightly overestimated 
within 3%, and the critical volume for a complete transition 
into the LS state is estimated to be about 53 A'^/f.u. 

the ambient structure, the energy preference of the HS 
state over the LS state is less than 0.15 eV/f.u. (see Fig. 
2 in Ref? ) , vifhereas the corresponding value we calculated 
is more than 0.5 eV/f.u. Therefore, the previous predic- 
tion of the HS-to-LS transition in BiCoOa with a small 
volume decrease of ~5% was overly optimistic.^ Instead, 
we find now that there is no readily switchable HS-LS 
transition is BiCoOa. 

As Bi^+ has a very similar ionic size as La'^^, it is rea- 
sonable to assume that when the local Co-0 bondlengths 
of BiCoOs in the HP PE phase become identical to those 
of LaCoOa in the LS state at low temperature (1.925 
A at 5 K)^ a complete transition into the LS state 
would be achieved in BiCoOa. This, together with the 
structural data of the 5.8 GPa phased allows us to esti- 
mate the critical volume to be about 52.8 A'^/f.u. (corre- 



sponding to a volume decrease of about 20%). Then, by 
extrapolating the eye-guided line of the V-P data points 
in the range of 2-6 GPa (see Fig. 1(b) in Refi), we may 
estimate the critical pressure to be about 8 GPa. It is 
important to note that the estimated critical volume is 
indeed well reproduced by our detailed calculations of 
the E-V curves (see Fig. 6, Vc ~ 53 A^/f.u.), which also 
nicely reproduce the equilibrium volume of the ambient 
phase within 3% and clearly indicate the HS-to-LS tran- 
sition via the mixed HS-I-LS state. This prediction of a 
complete transition into the LS state awaits a further HP 
study. 



IV. CONCLUSION 

To conclude, using configuration-state-constrained 
GGA-I-J7 calculations, we demonstrate that the giant 
tetragonal ferroelectric distortion of BiCoOa is driven by 
the strong Bi-0 covalency (rather than by lifting of the 
orbital degeneracy of the HS Co^+ ions) and is further 
stabilized by a subsequent xy-type 00. Moreover, our 
results show that the pressure induced HS-to-LS transi- 
tion is via a mixed HS-I-LS state, and that the transition 
would be complete upon a large volume decrease of about 
20% (under about 8 GPa). The mixed HS+LS state well 
accounts for the available experimental resultsi^ 
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